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1 Introduction 

> 

In these lectures I shall describe some of the latest data on particle production in high- 
ly-) . energy collisions and compare them with theoretical calculations and models based on 
QCD. The discussion will concentrate mainly on hadron distributions in jets, which are 
the manifestation at the hadronic level of hard (high-momentum-transfer) scattering of the 
Q\ ; partons (quarks and gluons) which are the fundamental fields of QCD. 
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In sect. |2|, the connection between parton and hadron distributions is made more precise 
using the concept of fragmentation functions. I concentrate in particular on the region of 
^| small momentum fractions, where interesting characteristic features of QCD are manifest. 

Next, in sect. |3], the various available models for the conversion of partons into hadrons are 
reviewed. In sect. |], the predictions of theory and models are compared with experimental 
data. After that, sect. |5] focuses on new data that show clearly the differences between jets 
that originate from quark and gluon fragmentation. 

Deep inelastic lepton scattering (DIS) at HERA is a copious source of jets; sect. |] dis- 
cusses new results on the properties of jets in the so-called current and target fragmentation 
regions. 

In sect. |7] 1 discuss new data on the fragmentation of heavy (b) quark jets into B mesons, 
and finally sect. |8] draws some brief conclusions. 

Many of the topics mentioned here are discussed more fully in ref. [[[]. In order to 
bring the discussion up to date, I have tried wherever possible to refer to the very latest 
experimental data. Therefore many of the references and figures concern preliminary data 
shown only at conferences, in particular at the International Europhysics Conference on 
High Energy Physics (EPS-HEP 99) held in Tampere, Finland, in July 1999. The cited 
contributed papers can be found through the EPS-HEP 99 web page || or, in the case of 
the large collaborations, through the collaboration pages |3j-[[l(J. 

1 Lectures at International Summer School on Particle Production Spanning MeV and TeV Energies, 
Nijmegen, The Netherlands, August 1999. 
2 Permanent address. 
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2 Jet fragmentation — theory 



We start with the basic factorization structure [1.1] of the single-particle inclusive distribu- 
tion, e.g. in e + e~ — > hX (fig. [I]): 

F h ( Xj s) = Y, -CfaaeWDiW*,*) , (!) 

■ J x Z 

i 

s = q 2 , x = 2p h ■ q/q 2 = 2E h /E cm (2) 

where Cj are the coefficient functions for this particular process (including all selection cuts 
etc.) and D\ is the universal fragmentation function (12[ for parton % — > hadron h. 

Although universal, fragmentation functions are factorization scheme dependent [[11 



If one tries to calculate them in perturbation theory, one encounters divergences associ- 
ated with the propagation of partons over long distances. In reality, however, partons are 
confined and cannot travel long distances. The perturbative divergences can be collected 
into overall factors that are replaced by non-perturbative factors taken from experiment. 
In this way the incorrect long-distance behaviour of perturbation theory is replaced by the 
correct long-distance features of QCD. However, the factorization of the divergent terms is 
ambiguous: one can choose to include different finite parts as well. This is the factorization 
scheme ambiguity. To specify the scheme requires calculation of the coefficient functions 
to (at least) next-to-leading order. This has only been done in a few cases. Thus there is 
need for theoretical work to make full use of the data on fragmentation functions. 

In certain kinematic regions, higher-order corrections are enhanced by large logarithms, 
which need to be resummed to all orders. Large logarithms of ratios of invariants may 
appear inside the coefficient functions Cj, for example in multi-jet processes when the 
angles between jets become small. In some cases these can be absorbed into a change of 
scale in the fragmentation functions. Examples of this will be encountered in sects. f| and 

The fragmentation functions D\ are not perturbatively calculable but their s-dependence 



(scaling violation) is given by the DGLAP evolution equation [|13|, [14], [15 

r) fl rj ' y 

Sy s D?(x, *) = £ j( i p d*, <*s(8))D}(x/z, s) (3) 

where Pji is the parton i — > j splitting function. Thus fragmentation functions can be 
parametrized at some fixed scale so and then predicted at other energies [12 . 

The most common strategy for solving the DGLAP equation is to take moments (Mellin 
transforms) with respect to x: 

D(N, s) = f 1 dx x N ^ D(x, s) , (4) 
Jo 
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the inverse Mellin transform being 

D(x,s) = — I dN x~ N D(N,s) , (5) 

where the integration contour C in the complex N plane is parallel to the imaginary axis 
and to the right of all singularities of the integrand. After Mellin transformation, the 
convolution on the right-hand side of eq. (^) becomes simply a product. 

The moments Pji of the splitting functions are called anomalous dimensions, usually 
denoted by jji(N,a s ). They have perturbative expansions of the form 



lAN,a s ) = f:^\N)(^-) n+1 . (6) 



We can consider fragmentation function combinations which are non-singlet in flavour, 
such as Dy = D qi — D q . or D qi — D qj . In these combinations the mixing with the flavour- 
singlet gluon drops out and for a fixed value of a s the solution is simply 

Dy(N,s)=Dy(N,S )(-) . (7) 

\Sq/ 

For a running coupling a s (s), the scaling violation is no longer power-behaved in s. The 
lowest-order form of the running coupling is 

as(s) = bwTW) (8) 

with b = (11C \ — 2n/)/127r, where Ca = 3 for QCD and n/ is the number of 'active' quark 
flavours, i.e. the number with m q -C s. Using this we find the solution 

Dy(N,s) = D v (N,s )(^j , dUN) = ^^, (9) 

which varies like a power of Ins. 

For the singlet fragmentation function 

£>5 = £(*>« (io) 

i 

we have mixing with the fragmentation of the gluon and the evolution equation becomes a 
matrix relation of the form 

d(Ds\_( 7« 2n /7OT \ / D s \ 

ds{D 9 )-{ lqg lgg ){ D g ) ■ ^ lL) 
The anomalous dimension matrix in this equation has two real eigenvalues 7± given by 

7± = | [igg + Iqq ± \l (igg ~ Iqqf + 8n ngqlqg ] . (12) 

Expressing Ds and D g as linear combinations of the corresponding eigenvectors D + and 
D_, we find that they evolve as superpositions of terms of the form (g) with 7+ and 7_ in 
the place of ^ qq . 
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At small x, corresponding to N — > 1, the q — > g and g — > (7 anomalous dimensions have 
a singularity, 

^-^T) + °( a ' ) - ^^^T) + °( Q S (13 > 

(Cf = 4/3), and we find 7+ — * 7 93 — > 00, 7_ — > ^ qq — > 0. Thus the low-x region requires 
special treatment, as we discuss in the following subsection. 



2.1 Small-x fragmentation 



At small x, multiple soft parton emission gives rise to terms enhanced by up to two powers 
of lnx for each power of a s . The leading enhanced terms can be resummed by changing 
the DGLAP equation (||) to 



d_ 

ds 



Di{x,s) = W 1 -P Jt (z,a s (s))D^(x/z,z 2 s) . 

■ J x Z 



(14) 



The fact that the scale on the right-hand side should be z 2 s rather than s follows from 
angular ordering of successive parton emissions [TJ], |Tj] . 



For simplicity, consider first the solution of eq. (|14]) for gluon fragmentation, taking a s 
fixed and neglecting the sum over different partons. Then taking moments as before we 
have 

d ~ 



dzz N ~ 1 P(z,a s )D(N,z 2 s) . 

Now if we try a solution of the form 

D(N, s) oc s ^ N ' as) 

we find that the anomalous dimension 7 (A/ - , a s ) must satisfy the implicit equation 

7 (JV, « s ) = f 1 dz z N-Wy(N,as)p( Z) ag ) . ( 17 ) 



(15) 
(16) 



When AT — 1 is not small, we can neglect the 27 (N, a s ) in the exponent of eq. (|17D and 
then we obtain the usual explicit formula for the anomalous dimension. For N ~ 1, the 
region we are interested in, the z — > behaviour P gg — > Caois/^z dominates, which implies 
that near N = 1 

lgg{N, a s ) = % S Ar , , I TTr I i (18) 



and hence 



7 99 (A^,a 6 



7T AT-l + 27 99 (AT,a s ) 



71 



\G A a 



1 



1 



2tt 



2tt 4^ 1 ' ) + 32VC , A a s 



(A^-l) 2 + 



(19) 



Thus for N 1 the gluon-gluon anomalous dimension behaves like the square root of 
a s . How can this behaviour emerge from perturbation theory, which deals only in integer 



A 



powers of a s ? The answer is that at any fixed N ^ 1 we can expand eq. fll^D in a different 
way for sufficiently small a s : 



lgg( N ,®s 



C A a s 



7T(N 



Ca® s 



tt J (N-iy 



+ ■ 



(20) 



This series displays the terms that are most singular as iV — > 1 in each order. These 
terms have been resummed in the expression (|T9"D, allowing the perturbation series to 
be analytically continued outside its circle of convergence \a s \ < (it/8Ca)\N — 1| 2 . By 
definition, the behaviour outside this circle (in particular, at AT = 1) cannot be represented 
by the series any more, even though it is fully implied by it. 

At sufficiently small x, the A" — > 1 singularity of the gluon-gluon anomalous dimension 
dominates in all fragmentation functions, and this in turn determines the asymptotic be- 
haviour of the single-particle inclusive distributions F h in eq. (|]). To predict this behaviour 
quantitatively we need to take account of the running of a s , which can be done by writing 
eq. (|16|) in the form 

/at ^ds r 

lgg\ N i "s) — 



D(N, s) ~ exp 



and noting that a s in the integrand should be a s (s'). We then use eq. (H) to write 



lgg{N, a £ 



W) d 4 = 4 



a s (s) 



Os 



da* 



at 



(21) 



(22) 



and hence 



D(N,s) oc exp - 




(N-l) 



+ 7Zl\br-AN - I) 2 + 



as=as{s) 



(23) 



The value of D(N,s) at A^ = 1 is simply the integral of the fragmentation function, 
which gives the average multiplicity, 



2C, 



2C A 
exp \ — — in s 
V no 



(24) 



Thus the average multiplicity of any hadron species should increase asymptotically faster 
than any power of Ins but slower than any positive power of s. Furthermore the relations 
([Uf) imply that the average multiplicities in gluon and quark jets are asymptotically in the 
ratio of their 'colour charges' C A and Cp [II? 



NMs) 



Ca 
C F 



9 
4 



(25) 



The behaviour of D(N,s) near N = 1 determines the form of small-x fragmentation 
functions via the inverse Mellin transformation (|5|). Keeping the first three terms in the 



Taylor expansion of the exponent, as displayed in eq. (p3|), gives a simple Gaussian function 
of N which transforms into a Gaussian in the variable £ = ln(l/x): 



xD(x, s) oc exp 



2a 2 



where the peak position is 



1 



,p 46a s (s) 4 
and the width of the distribution of £ is 



In s 



(26) 



(27) 



2^ 



(7 



246 V C A a£(s) 



oc (Ins 



(2J 



Thus the effect of resummation is to generate a characteristic hump-backed shape in 
the variable £ = ln(l/x), with a peak that moves up and expands slowly with increasing 
s. Including also next-to-leading logarithms, one obtains what is commonly known as the 
modified leading-logarithmic approximation (MLLA) Jh" 
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3 Hadronization Models 



3.1 General ideas 



Before discussing specific models for the hadronization process, we should review some 
general ideas that have proved useful in interpreting hadronization data. 

• Local parton-hadron duality p2 |. Hadronization is long-distance process, involving 



only small momentum transfers. Hence the flows of energy-momentum and flavour 
quantum numbers at hadron level should follow those at parton level. Results on 
inclusive spectra and multiplicities support this hypothesis. 



Universal low-scale a s [§3], [24], |25|. Perturbation theory works well down to low 
scales, Q ~ 1 GeV. Assume therefore that a s (Q 2 ) can be defined non-perturbatively 
for all Q, and use it in evaluation of Feynman graphs. This approach gives a good 
description of heavy quark spectra and event shapes. 



3.2 Specific models 

The above general ideas do not try to describe the mechanism of hadron formation. For this 
we must resort to models. The main current models are cluster and string hadronization. 
We describe briefly the versions used in the HERWIG and JETSET Monte Carlo event 
generators, respectively. In both parton shower initiated by the hard process evolves 

perturbatively, according to the DGLAP equation, until the scale of parton virtualities has 
fallen to some low value Qq ~ 1 GeV, whereupon the non-perturbative processes assumed 
in the model take over (fig. |2[). 



f, 



Figure 2: Cluster and string hadronization models. 



Cluster model [26|-[29|. The model starts with non-perturbative splitting of all gluons 
after the parton shower, g — > qq. Colour-singlet qq combinations have lower masses 



and a universal spectrum due to the preconfinement |3(], [Hp property of the shower 
(fig. [3] P2"l). These colour-singlet combinations are assumed to form clusters, which 
mostly undergo simple isotropic decay into pairs of hadrons, chosen according to the 
density of states with appropriate quantum numbers P7fl . This model has few param- 
eters and a natural mechanism for generating transverse momenta and suppressing 
heavy particle production in hadronization. However, it has problems in dealing with 
the decay of very massive clusters, and in adequately suppressing baryon and heavy 
quark production. 



String model |p3|- ||36|| . This model is based on the dynamics of a relativistic string, 
representing the colour flux stretched between the initial qq. The string produces a 
linear confinement potential and an area law for matrix elements: 



\M(qq-+ h---h n )\ 2 oce 



-bA 



(29) 



where A is the space-time area swept out (fig. 0). The string breaks up into hadrons 
via qq pair production in its intense colour field. Gluons produced in the parton 
shower give rise to 'kinks' on the string. The model has extra parameters for the 
transverse momentum distribution and heavy particle suppression. It has some prob- 
lems describing baryon production, but less than the cluster model. 



The UCLA model [57|, [38| is a variant of the JETSET string model which takes the 
above area law for matrix elements more seriously, using it to determine the relative 
rates of production of different hadron species. This results in heavy particle suppres- 
sion without extra parameters, the mass-squared of a hadron being proportional to 
its space-time area. At present the model still uses extra parameters for p? spectra, 
and again has some problems describing baryon production. 
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;ure 3: Cluster model: mass distribution of qq pairs. 




Figure 4: String model: space-time picture. 
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Table 1: Meson yields in Z° decay. Experiments: A=Aleph, D=Delphi, L=L3, M=Mark 
II, 0=Opal. Bold: new data this year. Underlined: disagreement with data by more than 



4 Single-particle yields and spectra 

Tables [1] and ^\ compare predictions of the above models^ with data on Z° decay from 
LEP and SLC. Of course, the models have tunable parameters, but the overall agreement 
is encouraging. As stated earlier, the main problems are in the baryon sector, especially 
for HERWIG. 

It is remarkable that most measured yields (except for the 0~ mesons, which have special 
status as Goldstone bosons) lie on the family of curves 

(n) = a(2J + l)e~ M/T (30) 



where M is the mass and T ~ 100 MeV (fig. | g(|). This suggests that mass, rather 



than quantum numbers, is the primary factor in determining production rates. Note that, 



3 Updated from ref. |§. 

4 Recent ALEPH HERWIG tuning with strangeness suppression 0.8 p9| 
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Table 2: Baryon yields in Z° decay. Legend as in table 1. 



surprisingly, the orbitally-excited J = | baryon A(1520) (not yet included in models) is 



produced almost as much as the unexcited J = | baryon E(1385) f|2 



At other energies, model predictions for identified particle yields are in broad agreement 
with e + e~ data (fig. | |]43|), but statistics are of course poorer. Charged particle spectra at 
low x agree well with the resummed (MLLA) predictions discussed in sect. |] over a wide 
energy range, as illustrated in fig. |7| 

In pp — > dijets the relevant scale is taken to be Q = Mjj sin 9 where Mjj is the 
dijet mass and 9 is the jet cone angle (fig. |8|). Results are then in striking agreement with 
MLLA predictions and with data from e + e~ annihilation at Q = ^/s (fig. |9|). 

New SLD data include hadron spectra in light quark (rather than antiquark) fragmen- 
tation, selected by hemisphere using the SLC beam polarization |4(J. One sees strong 
particle/antiparticle differences in the expected directions (fig. [TUp, bearing in mind the 
predominance of down-type quarks in Z° decay. 



5 Quark and gluon jets 



DELPHI [47j select gluon jets by anti-tagging heavy quark jets in 'Y' and 'Mercedes' 
three-jet events (fig. [11]) . As expected, the higher colour charge of the gluon (Ca = 3 vs. 
Cp = 4/3) leads to a softer spectrum and higher overall multiplicity (fig. [12]). In general 
the relative multiplicities of identified particles are consistent with those of all charged, 
with no clear excess of any species in gluon jets (fig. [13]). In particular there is no enhanced 
0(1020) or r] production: 

DELPHI 0: N g ((j))/Ng((j)) = 0.7 ± 0.3 
OPAL @: N g (r})/N q (r}) = 1.29 ±0.11 



OPAL [j|9| select gluon jets recoiling against two tagged 6-jets in the same hemisphere. 
Monte Carlo studies indicate that such jets should be similar to those emitted by a point 
source of gluon pairs. The qualitative message from the data is again clear (fig. |T^): Gluon 
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Figure 6: Particle yields in e + e annihilation. 
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Figure 7: Low-x fragmentation in e + e annihilation. 




Figure 8: Cone angle in pp — > dijets. 
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Figure 9: Low-x fragmentation in pp — > dijets. 
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Figure 10: Normalized particle-antiparticle differences in quark jet fragmentation. 
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Figure 11: Selection of gluon jets by DELPHI. 




Figure 12: Charged particle spectra in quark and gluon jets. 
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Figure 13: Comparisons of particle spectra in quark and gluon jets. 
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Figure 14: Momentum fraction and rapidity distributions in quark and gluon jets. 



jets have softer fragmentation than light quark jets, and higher multiplicity. The precision 
of the data is now such that next-to-leading order calculations of the relevant coefficient 
functions, taking into account the experimental selection procedures, are needed to check 
universality of the extracted gluon fragmentation function. 

The ratio of gluon/quark multiplicities at low rapidity (large angle) is close to the ratio 
of colour charges r = Ca/Cf = 2.25, in agreement with local parton-hadron duality: 

OPAL: r ch {\y\ < 1) = 1.919 ± 0.047 ± 0.095 . 

According to eq. (0), the overall multiplicity ratio should also approach the value Ca/Cf 
asymptotically, but at present energies the contribution from higher rapidities is substantial 
and this leads to a smaller ratio. 



Monte Carlo studies [||J suggest that a better measure of Ca/ Cp is obtained by selecting 
low-momentum hadrons with relatively large transverse momentum (i.e. low rapidity). This 
gives 

OPAL: r ch (p < 4, 0.8 < p T < 3 GeV) = 2.29 ± 0.09 ± 0.015 . 

DELPHI [|5(J have observed scaling violation in quark and gluon jet fragmentation 
separately (fig. |1^) by studying the dependence on the scale 

k h = sin(0/2) ~ \^yl (31) 



where 9 is the angle to the closest jet and y% is the Durham jet resolution [51] at which 3 



jets are just resolved. This is expected to be the relevant scale when y% becomes small. One 



sees clearly that there is more scaling violation in gluon jets (fig. (T6J) . The ratio provides 
another measure of Ca/Cf- 

DELPHI: r scMo i = 2.23 ± 0.09 ± 0.06 . 

A crucial point in the DELPHI analysis is that 3-jet events are not selected using a 
fixed jet resolution y cnt , but rather each event is clustered to precisely 3 jets. This avoids 
'biasing' the gluon jet sample by preventing further jet emission above y cut . 
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Figure 15: Scale dependence of quark and gluon fragmentation. 




Figure 16: Logarithmic gradients of quark and gluon fragmentation. 
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Figure 17: Rapidity plateau in two-jet events. 




Figure 18: Average qq and gg multiplicities deduced from 2- and 3-jet events. 



The same point is well illustrated in analyses of average multiplicities in 2- and 3-jet 
events |5^, 54|. If N qq -(s) is the 'unbiased' qq multiplicity, then in events with precisely 
2 jets at resolution y CVLt there is a rapidity plateau of length ln(l/y cut ) (see fig. [LTD and the 
multiplicity is 

N 2 {s, y CVLt ) ~ N gq {s ycat ) + ln(l/y cut )A^(sy cut ) (32) 

where N'(s) = sdN/ds. Clustering each event to 3 jets we get this multiplicity with y 3 in 
place of y cut , plus an unbiased gluon jet: 



N 3 (s) ~ N 2 (s,y 3 ) + ±N gg (sy 3 ) . 



(33) 



Thus one can extract the unbiased gg multiplicity, plotted in fig. |T8 
The ratio of gg/qq slopes gives yet another measure of Ca/Cf |5 



vs. Pi ~ 



r mu i t = 2.246 ± 0.062(stat.) ± 0.080(sys.) ± 0.095(theo.) . 



6 Current and target fragmentation in DIS 



The HI and ZEUS [[57]] experiments at HERA have studied the distributions of x p = 
2\p\/Q in the current and target hemispheres in the Breit frame of reference. The Breit 



1 v 



p, 




(l-x)Q/2x 



CURRENT TARGET 
Figure 19: Breit frame current and target regions in DIS. 
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Figure 20: Contributions to the final state in DIS. 



frame is the one in which the 4-momentum of the virtual photon exchanged in DIS lies 
entirely along the negative z-axis, = (0,0,0, —Q), while the target proton has = 
(Q, 0, 0, Q)/2x (neglecting the proton mass), where x = —q 2 /2p ■ q is the Bjorken variable. 
In this frame, to zeroth order in a s , the virtual photon simply strikes a constituent of the 
target with momentum Q/2 and reverses its momentum. The remnant of the target is then 
left with momentum (1 — x)Q/2x (fig. [T9| ) . 

In higher orders one expects the current hemisphere to contain fragmentation products 
of the current jet (C in fig. |2(|), similar to half an e + e~ event. In the target hemisphere, 



the contribution Tl is similar to C, T2 gives extra particles with x p < 1, while T3 gives 
Xpkjl, generally outside detector acceptance. 



In the current hemisphere the charged multiplicity is indeed similar to e + e (fig. |21 



fjTjj). Differences at low Q 2 are consistent with the expected boson-gluon fusion 



contribution. The distribution of £ = ln(l/x p ) is also similar to e + e , i.e. close to 
Gaussian with little Bjorken x dependence (fig. |22"|) . 

At low Q 2 there is evidence of strong subleading corrections. The distribution is 
skewed towards higher values of £ (smaller x p ), contrary to MLLA predictions (fig. 



1 Q 
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Figure 21: Charged multiplicity in current hemisphere. 
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Figure 22: Fragmentation in DIS. Upper data (heavy curve) target region, lower data (light 
curve) current region. 
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Figure 23: Skewness in current fragmentation region. 



The quantity plotted is 



Skewness = /(£ - £) 3 > / ((£ - 



3 

2\2 



(34) 



On the other hand, the data lie well below the fixed-order perturbative prediction 
at low x p and Q 2 (fig. El). Discrepancies could be due to power-suppressed (l/Q 2 ) 



corrections, of dynamical and/or kinematical origin. The bands in fig. 
to an ad-hoc correction factor 



M\ correspond 



(0.1 < m cS < 1 GeV). 



(35) 



In the target hemisphere there is also disagreement with MLLA J57|, possibly due 
to the T3 contribution "leaking" into the region x p < 1. If anything, Monte Carlo 
models predict too much leakage (fig. pE|). Little Q 2 dependence is evident. 



7 Heavy quark fragmentation 



New data on b — > B fragmentation from SLD p5 |, using high-precision vertexing, discrim- 
inate between parton-shower plus hadronization models (fig. p6|). Note that the data have 
not yet been corrected for detector effects. 

In general one expects the b quark to lose energy by gluon emission in the parton shower, 
and then to suffer a further non-perturbative energy loss during hadronization. The latter 
is conventionally parametrized by convoluting with the Peterson function ||60|| : 

/(*) = -(l--- -^-) 2 (z = x B /x b ) (36) 



z \ z 



on 




Figure 24: Scaling violation in DIS fragmentation. 
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gure 25: Target fragmentation compared with models. 




Figure 26: SLD data on b — > B fragmentation compared with models. 



Including more perturbative QCD leads to a reduction in the amount of non-perturbative 
smearing required to fit the data, and hence to a smaller fitted value of the Peterson 
parameter e&: 

Pure Peterson | 59fl : e& = 0.036. 

0.006. 



JETSET parton shower + Peterson 



0.002. 



NLLA QCD + Peterson [§T]: e b 

Here NLLA (next-to-leading logarithmic approximation) refers to an analytical pertur- 
bative calculation that goes beyond the parton shower approximation. The calculation of 
ref. [|61 



was fitted to earlier ALEPH data |p2| , as shown in fig. |27. 

In the universal low-scale a s model, the perturbative prediction is extrapolated smoothly 
to the non-perturbative region, with no Peterson function at all |24[ . 



8 Conclusions 

We have seen that experimental studies of fragmentation are yielding large amounts of 
new data for comparison with theoretical predictions and models. Especially impressive 
is the success of the MLLA perturbative predictions in accounting for the general shape 
and energy-dependence of fragmentation at small momentum fractions. We are now at the 
stage when more detailed comparisons await new theoretical input, in the form of coefficient 
functions that take account of selection procedures, especially for gluon jets in e + e~ final 
states. 

Comparisons between data and hadronization models suggest that particle masses, 
rather than quantum numbers, are the dominant factor in suppressing heavy particle pro- 
duction. Baryon production is not yet well described by any model. 





Figure 27: ALEPH data on b — > B fragmentation compared with NLLA QCD. 



Quark and gluon jets have the expected differences and these can be used to measure 
the ratio of colour factors Ca/Cf- There is no strong evidence yet for different particle 
content in gluon jets. 

Fragmentation in deep inelastic lepton scattering shows some disagreements with per- 
turbative predictions. It is not yet clear whether these are due to higher-order or non- 
perturbative effects. 

New precise b quark fragmentation data from Z° decay are now available and put strong 
constraints on models for heavy quark hadronization. The data suggest that perturbative 
effects dominate. 
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